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ارتقای کارابی شناسایی پارامترهای فیزیکی سازه با استفاده از کاهش تکینگی ماتریس‌های پاسخ در حل معکوس 
معادلات ح رکت * 


(۳) 5 0) 


مجید قاسمی بابک تقوی"" سعید عباس بندی 
چکیده در اين مقاله شناسایی ماتریس‌های مشخصه سیستم حطی (شامل ماتریس‌های جرم. میرایی و سختی) بااستفاده از حل معکوس 
معادلات حرکت, با تم رکز بر کاهش شرایط بدوضعی دستگاه معادلات حاصل,ء در حوزهُ فرکانس ارائه گردیده است. این روش با اندازه‌کیری 
پاسخ‌های شتاب در درجات آزادی سازه تحت یک تحریک اجباری, به شناسایی سیستم و تشخیص خرابی‌های محتمل می‌پردازد. بدوضعی 
معادلات حرکت و وجود نویز در پاسخ‌های شتاب اندازه‌گیری‌شده به‌عنوان یک مشکل ناگزیر» موجب ناپایداری پاسخ سل تشخحیص 
سیستم, خحطای فزاینده و عدم قابلیت اعتماد در ماتریس‌های مشخصه تعیین‌شده می‌شود. دراین مقاله الکوریتمی برای بهبود بحشیدن به مشکل 
بدوضعی حل معکوس مسئله در روند شناسایی سیستم ارائه گردیده است که یک روش خحاص بالامثلثی‌سازی ماتریس است. الکوریتم 
پیشنهادی با شناسایی بردارهای موازی و شبه‌موازی موجود در ماتریس ضرایب دستگاه معادلات خعطی به حذف بردارهای وابستة حطی 
می‌پردازد و با کاهش تکینکی ماتریس. موجب پایدارسازی و ارتقای پاسخ می‌گردد. در تخمین بهینة ماتریس‌های مشخحصه سیستم در روند 
شناسایی از روش حداقل مربعات و تابع جریمه استفاده شده است. به‌منظور بررسی کارایی الگوریتم پیشنهادی, ماتریس‌های مشخحصه قاب 
هشت طبفهٌ غیربرشی با میرایی نامتناسب با استفاده از الگوریتم کاهش تکینگی در شناسایی سیستم‌تعیین شده‌اند. نتایج نشان می‌دهد استفاده 
از الگوریتم پیشنهادی در بهبود تخمین و پایداری پاسخ‌ها کاملا مفید است. 


واژه های کلیدی شناسایی سیستم سازه. حل معکوس معادلةً حرکت. مسئله بدوضع. تکینگی ماتریس. 
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* تاریخ دریافت مقاله ۹۳/4/۳۱ و تاریخ پذیرش آن ۹۶/4/۲6 می‌باشد. 
(۱) استادیان دانشگاه آزاد اسلامی. واحد قزوین, دانشکده عمران و نقشه‌برداری» گروه عمران. قزوین. 


(۲) نویسنده مسئول: دانش آموخته کارشناسی ارشد دانشگاه آزاد اسلامی» واحد قزوین. دانشكدة عمران و نقشه‌برداری» گروه عمران قزوین. 
ومع 2027108].ظ 
(۳) استاد. دانشگاه بین‌المللی امام خمینی(ره؛ دانشکدة علوم پایه. گروه ریاضی کاربردی» قزوین. 
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معد مه 
امروزه آزمایش غیرمخرب سازه‌ها (۳1) با هدف 
شناسایی سیستم. تشخیص خرابی و پایش سلامت 
سازه در بسیاری از سازه‌های مهم مورد استفاده قرار 
می‌گیرد [1]. اين روش‌ها عموماً بر پایة الگوریتم‌های 
شناسایی سیستم استوار هستند و در آنها به‌کمک 
اندازه گیری‌های میدانی. ساختار و پارامترهای مدل 
دینامیکی شناسایی می‌شوند. روند شناسایی سیستم 
مبتنی بر حل معکوس معادلات حرکت است که در آن 
پارامترهای سیستم به کمک اندازه‌گیری سیگنال‌های 


ورودی و خروجی آميخته با نویز تعیین می‌گردد. 

به‌طور کلی مسائل معکوس را می‌توان جزء 
بارزترین مسائل فیزیکی دانست که به لحاظ ریاضی 
بدوضم هستند مسائلی مانند مسائل معکوس پراکنش 
دزد کمن شنک اقا هه از ادها کی انز 
فاز امواج پراکنده‌شده به‌وسیلةٌ شیء تصویربرداری 
به‌روش توموگرافی در مبحث پزشکی. مسائل معکوس 
استیفان برای مدل‌سازی ذوب یخچال‌های قطبی. 
نمونه‌هایی از مسائل فیزیکی بدوضع هستند که مبتنی 
بر حل معکوس می‌باشند [2]. 

در حبل عددی دستگاه‌های معادلات. بنابر 
وضعیت جواب مسئله یعنی وجود. یکتایی و یا 
وضعیت پایداری جواب. روش‌های گوناگونی مورد 
استفاده قرار هی کیبل به‌عبارت دیگخا زمانی که سه 
قرط زا زمره( کت ای( پایدار جزات 
تأمین باشند روش‌های معمول حل دستگاه‌های 
معادلات می‌توانند مورد استفاده قرار گیرند اما در 
صورت عدم برقراری هر یک از شروط یاد شده با یک 
مسئلٌ بدوضع روبه‌رو خواهیم بود که حل آن نیازمند 
ردان و خوه ات در اران رن 6 خاداساره 
((1865-1963) ۳۲۱۵0۵۴2۲0 ومصصماقو ععا0ع12) او لین 
کسی بود که بحث دستگاه‌های معادلات خوش‌وضع و 
بدوضع را مطرح نمود. وی دستگاه‌های معادلاتی را که 
در آنهاسه شرط یاد شده برقرارند را مسائل 


نشریهُ مهندسی عمران فردوسی 


تا کارا ای ماه قزر یکین شاری ۲ شاه از 


خوش‌وضع (صعاطامتم صمتاتع0) [[۳۵)» و مسائلی را 
که در آنها یکی از سه شرط نقض گردیده باشد. مسائل 
«بدوضع(صعادمیج همناندی )» نام گذاری کتر ق 
البته هادامارد بر این عقیده بود که مسائل بدوضع تنها 
و ی تمس فساقام ق کی 
موضوعیت یابند[2]. به‌طور کلی مسائل بدوضع به دو 
دسته طبقه‌بندی می‌شوند. مسائل دارای «نقص رتبه» که 
در آنهاماتریس ضرایب دستگاه معادلات دارای 
تعدادی مقادیر تکین کوچک است به‌طوری که مقدار 
آنها بسیار کوچک‌تر از سایر مقادیر تکین این ماتریس 
است. نوع دیگر از مسائل بدوضع مربوط به حالتی 
است که رتبه و مرتبهٌ عددی ماتریس ضرایب دستگاه 
معادلات پر است. اما باز هم تفاوت آشکاری بین 
مقادیر تکین کوچک و بزرگ این ماتریس وجود دارد 
که این حالت تحت عنوان «مسائل بدحالت ( 00560 ۲11 
حع[ا0۲۵)) شناخته می‌شود [3]. بر این اساس رامپ 
(۳۳۵ 5.۷۲) [4,5] روشی برای معکوس کردن 
ماتریس بدوضع با عدد وضعیت دل‌خواه ابداع نمود. 
ویژگی روش او این بود که محاسبات با دقت دلخواه 
فقط در بخشی از کل فرآیند حل مورد نیاز بود. کارایی 
روش رامپ در مقالات دیگر [6-8] نیز بررسی شد. 
اوهتا و همکاران ([2 6 012 .1) [9] روش رامپ را با 
کمی تغییرات در تعداد ارقام با معنی محاسباتی بررسی 
کردند. تعداد زیادی از روش‌های پيشنهادشده برای حل 
دستگاه‌های بدوضع به لحاظ تعداد فرآیند تکراری 
مورد نیاز برای هم‌گرایی مقرون‌به‌صرفه نیستند 
(پرهزینه هستند) اگرچه برخحی از آنها مانند روش 
مورینسکی (۷]0527051 ,>) [10] و اسستراند 
0.50007 [11] از دقت مناسبی نیز برخوردارند. 
روش دیگری که برای کاهش مشکل بدوضعی پیشنهاد 
شده است الگوریتمی موازی برای تبدیل هاوسهولدر 
است که توسط لو و همکاران (21 6۶ دارآ .۷۲ [12] 
توسعه داده شده است. استفاده از تجزیه مقادیر 
تکین (8۷) در حل دستگاه معادلات بدوضع اطتین 


سال بیست و هشتم شمارهٌ دو ۱۳۹۵ 


مجید قاسمی- بابک تفوی-سعید عیاس بندی 


یک روند طولانی توسط چند ریاضیدان برجسته توسعه 
پات در ایا تروش اما ری فان ویتکا 
معادلات به بردارهای ویژه و مقادیر ویژهُ خود تجزیه 
فین گردنت. تسیز معسادلات :به‌دسستت: آمداه بر آی مل 
دستگاه مورد استفاده قرار ,کی تنل اگرچه ٩۷۲‏ 
به‌عنوان مهم‌ترین روش پایدارساز مسائل بدوضح 
شتاحته می‌شود. بااین حال دارای اشکالاتی یز 
می‌باشد مانند از دست دادن برخی اطلاعات مفید در 
همان گام اول این روش (به‌دست آوردن مقادیر تکین 
تاو ال ان تفا اسان ورگ متا 
بسیاری ارائه نموده‌اند که از آن جمله می‌توان به 
الگوریتم دوقطری‌ساز گولاب و کاهان ( للع 614 
م1 220) [13] الگوریتم بوزینگر و گولاب ( ۳.۸ 
طبامی 06.۳ 20 عومزون:ظ۳) [14] الگوریتم گولاب 


کاهان_رینسک (0.۳۴۵10600) 0ص حاتتامی 6:11) [15] 
اشاره نمود. الگوریتم‌های دیگری نیز بر پایُ روش‌های 
ذکرشده توسط ۳ جان (0۳27 1.۳) (۷1-جان) 
[16] دمل و کاهان (1۵12۳0 ۲۷۰ 290 [عصجع(1 ل) [17] 
و فرناندو و پارلت (۴۵:۱66 5.۲ 200 ۳6۳8۵800 1>,۷) 
(الگوریتم0۲05) [18] ارائه گردید. هم‌چنین کاربرد 
0 در حل برخی دستگاه‌هایی که دارای مقادیر یه 
بسیار کوچک هستند دچار مشکل می‌شود. برای حل 
تن مشکل هنسن (۴120568 6) [19] و واراه ( 3.1۷۲ 
[20,21] روش تجزيء ناقص مقادیر تکین 
(151)را پيشتهاد بر دند: با توه به آهست سل 
قطع در روش تجزیه ناقص مقادیر تکین کوشش‌های 
زیادی برای پيشنهاد محل قطع بهینه انجام شده است. 
نمونه‌ای از این روش‌ها که کاربردهای زیادی نیز دارند 
روش 01۲۳۷۲6 و روش ۲095-7۵11021107 است که 
هردو توسط هنسن [22,23] توسعه یافته‌اند. ولوخ و 
و پلنی («مصاز ۷ 0 0جج ما۷۵ 06:۲۷ [24] نقطه قطع 
بهینه را برمبنای تصوير کردن مقادیر ویژهُ بسیار کوچک 
به فضای خالی ماتریس بردارهای ویژه محاسبه نمودند 


آنها نام روش خود را بهبوديافتةٌ تجزیه ناقص مقادیر 


سال بیست و هشتم شمارة دو ۱۳۹۵ 


۳۵ 


تکین (115۷۳) نامیدند. یکی دیگر از معمول‌ترین 
روش‌های حل مسائل بدوضع که یک روش منظم‌ساز 
شناخته می‌شود توسط تیخونوف (1110000۷ .ه) 
[25] ارائه گردید اين روش خصوصا در حل مسائل 
معکوس کاربرد گسترده‌ای دارد و در آن به‌منظور 
پایدارسازی مسئله اطلاعاتی در خصوص جواب مسئله 
تن اضافه می‌گردد. 

درحصوص استفاده از الگوریتم‌های پایدارساز 
مطرح‌شده در مسائل شناسایی سیستم سازه‌ای نیز 
تحقیقات فراوانی صورت گرفته است. لسلر و همکاران 
راو ]6 1,601676 26[604] شناسایی نیروهای موجود بر 
روی سازه را از طریق اندازه‌گیری‌های شتاب بررسی 
کردند. آنها در روند شناسایی روش خود از تجزیه 
ناقص مقادیر تکین با نقطه‌قطع تعیین‌شده در روش 
۲6 استفاده نمودند. در ادامه نردب رگ و 
گوستافسون (061815507 [) [27] در روند این 
شناسایی از منظم‌سازی تیخونوف استفاده کردند. دیکن 
و همکاران (21 61 1060 28[6۷] از روش منظم‌سازی 
تیخونوف برای شناسایی بردارهای ثاببت دینامیکی در 
نسازه‌های ها مرانی غیر خطتی, انسفاده نمودتاه نتو و 
همکارانش (64۵1 ۷ ۳) [29] با این ایده که حل 
مشکل بدوضعی و نقص رتبه مهم‌ترین مسئله در روند 
شناسایی پارامترهاست روشی منظم‌ساز برای به‌روز 
نم دن فلال تاره ارانه تمودنتت کته هو ان یش‌کادتن 
عددی ذکرشده در فرآیند بهینه‌سازی کاهش یافته بود. 
روش آنها براساس حداقل مربعات متغیرهای محدود 
(2۲0او و ]وهع,] ۷۵۲1۵۳0169 0۲۳۴060ظ) 5ب ۷ظ بود و 
در آن بر روی راهکارهای ریاضی حل عددی مسئله 
بهینه‌سازی تمرکز شده بود. آنها روش را بر روی یک 
سازه با ورودی و خروجی مشخص اعمال نمودند و در 
آن تابع پاسخ فرکانسی را به‌عنوان داده‌های مسئله 
انتخاب گردند. آنان نشان دادند روش ارائه‌شده دارای 
سرعت هم‌گرایی بالاتر و دقت قابل قبول بوده است. 
تیتیوروس و فریسول (۳۲۱5۷۷76۱۱ ]1۷۲ 200 دناب11 ظ) 
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[30] روشی را برای به‌روز نمودن مدل مبتنی بر 
حساسیت ارائه نمودند که در آن روش بهینه‌سازی 
براساس تجزیه مقادیر تکین بود. آنها بیشتر در مورد 
حل مسئله بهینه‌سازی و روش‌های متفاوت آن تحقیق 
ره و مه نی فا میگ تشه 
مختلفی را به‌کار بردند و نحوه هم‌گرایی محاسبات را 
در این توابع بررسی کردند. در انتها روش بر روی 
مد ها عددعم اعمال نی فان مورد بابک 
قرار گرفت. 

در اين مقاله یک الگوریتم جدید تحت عنوان 
«کاهش تکینگی ماتریس» برای حل مشکل بد وضعی 
مسله شتامایی میستم سازه‌ای ارانه گردیه اسست که 
از تکنیک‌های جبر خطی بهره می‌گیرد. روش شناسایی 
مورد استفاده در این مقاله با استفادهٌ مستقیم از معادلات 
حرکت و به‌کارگیری بهینه‌ساز تابع جریمه به شناسایی 
پارامترهای سیستم در حوزه فرکانس می‌پردازد. باید 
توجه داشت در عمل یک هم‌پوشانی بین روش‌های 
شاوی و سس مساو وعضر 1۳1۵33۵ 
الگوریتم ارانه‌شده در این مقاله ضمن کاهش 
تکینگی‌های ماتریس ضرایب باعث تأثیر بهتر بهینه‌ساز 
تابع جریمه می‌شود که به شکل یک پیش‌شرطساز در 
فرم دستگاه معادلات جبری نمایان گردیده است. در 
ی هیاس سا کی کی 
ماتریس» ضمن کاهش عدد وضعیت ماتریس ضرایب 
دستگاه معادلات دقت و پایداری پاسخ مسئلةٌ شناسایی 
افزایش یابد. کارایی و قابلیت الگوریتم مذکور بر روی 
دو سیستم سازه‌ای با تعداد طبقات متفاوت و 
بارگذاری‌های مختلف (هارمونیک (1020 ۲12۳00816). 
تصادفی (1020 حصمصعظ) و ظ۳۳ (ر 80عباوظ 
0 569۲6۸66 #تفصنظ صمل‌ه۳2)) بررسی ۳ دیده 
است. در این‌جا این نتایج برای یکی از اين سازه‌ها که 
دک شوه هرت اه پاام رای یانب ای ارات 


گردیده أبیت: 
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زا کارا سای پراش فا فیریکین سار کم اهاز 


روش مستفیم شناسایی سیستم در حوزه فرکانس 
و استفاده از ایدهٌ کاهش تکینگی ماتریس 
در روش مستقیم شناسایی» معادلهٌ تعادل حرکت سازه 
خطی با 9 درجه آزادی و فرض میرایی ویسکوز و با 
درنظر گرفتن هه در سضانیه از ادی از و گام‌های 
زمانی اندازه‌گیری شده (7 گام زمانی) منجر به 10 < 8 
معادله خواهد شد. با جداسازی ماتریس‌های 
خحصوصیات سیستم و پاسخ‌های سازه به‌عنوان 
مجهولات و معلومات شناسایی سیستم. فرم بلوکی 


1 
)۱( مرو < تلا اک مها ممرط]۷] 
1 


با تبدیل فوریه طرفین معادله و بیان پاسخ‌ها و 
نیروی ورودی در حوزه فرکانس معادله قابل بازنویسی 
است: 

جمدو(ه) لا 
دو(ه) ٩‏ < مرو( تا | جک مج مدو]۷] 
«مدو(صه) 1 

)۲( 
که (ه) (ه) (م)لا و (ه) ب‌ترتیب 
پاسخ‌های شتاب. سرعت. جابه‌جایی و نیروی ورودی 
در حوزهة فرک‌انس بوده و همگی مقادیر مختلط 
می‌باشند. در صورتی‌که مقادیر پاسخ‌ها در حوزهء زمان 
نویزدار باشند. مقادیر پاسخ‌های فرکانسی نیز نویزدار 
خواهند بود. ولی در موارد عملی معمولاً توزیع نویزها 
در همه ناحیه‌های فرکانسی یکسان نیست و میزان 
نسبت نویز به پاسخ فرکانسی در مقادیر فرکانس‌های 
خیلی کم و زیاد و هم‌چنین در فرکانس‌های 
آنتی‌رزنانس‌ها بیشتر از سایر نواحی خواهد بود. در این 
مقاله به‌منظور اجتناب از وابستگی روش نسبت به نحوة 
توزیع نویز در محیط فرکانس نویزها در حوزه زمان با 
توزیع نرمال ایجاد شده و به پاسخ‌ها اضافه شده‌اند و 


سپس پاسخ‌های نویزدار با تبدیل فوریه به پاسخ‌های 


سال بیست و هشتم شمارهٌ دو ۱۳۹۵ 
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قالش شلیا بان 

پا معرفسی مس اتریس 
اس لا لا مسسه تا 5 ۳ به‌عنوان 
۱ 
و "امک مسا رما < وروولا بسهعنسوان 
ب‌صورت زیر نوشته می‌شود که پاسخ آن جواب مسئلة 
شاشاتی فتاه 


۳( ۳ < ولا موس 


در معادلهٌ جبری فوق ماتریس یک ماتریس 
فترامعین با عنده وضعیت بسیار پالامسته بنابراین 
درصورتی که پاسخ سازه به نویز آغشته باشد» 
جواب‌های محاسبه‌شده برای 0) با استفاده از رابطه (۶) 
منجر به خطای بسیار زیادی خواهد شد. دراین‌صورت 
ماتریس نیروی ماندگار سازه در محیط فرکانس» ول 
در 7 درجه آزادی سازه و 7 گام زمانی به‌صورت زیر 
قابل بیان است: 
(( کر 

باتوجه به مختلط بودن مقادیر مت در حسل 


معادلات می‌توان از کمینه بودن مقادیر حقیقی و با 


مقادیر قسمت موهومی و به‌صورت حداگانه استفاده 


نمود. 
از 7 -_ ۳۳ 
(۵-لف) ۳( ال و3۳8 مب 
1 - 
(۵-ب) ت ۳ > ۳ ۳ 


که 0 و 4 قسمت‌های حقیقی و موهومی 
ماتریس‌های پاسخ و ۲ و ,۲ قسمت‌های حقیقی و 
موهومی ماتریس‌های نیروی ورودی در حوزه فرکانس 
هت غفوری آشئبانن و قامنمی [31] ناست کردتل 
می‌توان تخمین ماتریس »0 را به‌عنوان تخمین درستی 
اما یبن خرن اسان سم خرن رکه 


سال بیست و هشتم شمارة دو ۱۳۹۵ 


۳۷ 


در این روش بهینه‌ترین خحصوصیات سازه به‌گونه‌ای 
انتخاب می‌گردد که مجموع مربعات نیروه‌ای 
ماند کاردشر هه درخانت اراد و خر هه گام‌های 
زمانی منتخب کمینه گردد. در حالت بهینه‌سازی نامقید. 
جواب بهینه‌سازی به شرح زیر خواهد بود: 


ت05 


0( ۶ 2 ۸۵ < 20 ود 
1 


در حالتی که قیدهای تقارن ماتریس‌های میرایی و 
سختی و قطری بودن ماتریس جرم اعمال گردند» 
می‌توان از روش‌هایی مانند روش تابع جریمه در روند 


بهینه‌سازی استفاده نمود. در این حالت قیدهای مزبور 


به شکل زیر قابل بیان است: 
ژد 2۳01 ظ ک 1 10۲ 0 < 0 < زره 
0 > 1 > ۲1+ 10۲ مزمبز0 < زره 


]0۲ 2۳41 1 > 1 > 0 
(۷) 


در این حالت تابع هدف بهینه‌سازی به‌صورت زیر 


,2-رصمبز0 < زره 


قابل بازنویسی است: 
20 1 ۳ 
ی 1 2۵ هد پر "1 


+ [-ز-ز ا-ز 1-۱ 


۱ (ووبز ود - 0 


)۸( 
در رابطه فوق م1 به‌عنوان ضریب تابع جریمه 
بیانگر حد بزرگی تابع جریمه می‌باشد و مقدار آن 
پراساس تابع هدف تعیین تین کت و33 انتخاب مقادیر 
کوحک برای 8 باعت تایر اندک قیدها در روند 
بهینه‌سازی می‌گردد. از طرف دیگر انتخاب مقادیر 
بسیار بزرگ برای آن باعث اختلال در روند بهینه سازی 
خواهد شد. بنابراین در این مقاله مقدار «8 برابر با 
مقدار حداکثر ۳ در تابع هدف درنظر گرفته شده 
است که با ایجاد خطا در برقراری قیدها؛ مقدار 
جریمه‌ای با مقیاس تابع هدف اولیه به تابع جریمه 

افزوده گردد. 
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درنهایت در حالت بهینه‌سازی مقید جواب بهينة 
: ی ۹۳ 
تابع هدف با فرآیند مشتق‌گیری 0 - رو قابل تعیین 
0 


)0 ۱۳۵۰۸ - رود 


در رابطة )٩(‏ تابع جریمه به فرم پیش شرط‌ساز 9 
در معادلهةٌ جبری (۶) نمایان شده است. 


الگوریتم کاهش تکینگی ماتریس #ط»ت0ع) 
((۲۹۲۷۲) امه ۵ زاتتهاناوطاه 
ایدةُ اصلی این مقاله بر پایةٌ کاهش تکینگی ماتریس "3۳ 
استوار است که در آن از نوع خاصی از بالامثلئی‌سازی 
استفاده شده است. الگوریتم ارائه‌شده قادر است با 
شناسایی ستون‌های وابستة عددی و وابستهٌ خطی 
ماتریس "7 این ستون‌ها را از ماتریس حذف کند. 
درحقیقت هر بردار ستونی در دستگاه معادلات خحطی 
(۶) یک معادله است که در سطری از ماتریس 5 قرار 
دارد. در روند حل این دتگاه تعداد این معادلات 
بسیار زیاد است بنابراین بهینه‌ساز تابع جریمه به معادله 
اضافه شده است. تعبیر جبرخطی استفاده از تابع 
جریمه به فرم یک پیش‌شرطساز در معادل (4) مشاهده 
می‌گردد. از نگاه جبر خطی عددی عملکرد یک 
پیش‌شرط ساز در دستگاه معادلات خطی دوران ماتریس 
ضرایب است. به‌طور مثال در فرم دومعادله دومجهول 
(مطابق شکل تن عملکرد به نبیین نقطة قطع 
وی کین بایان تن پاینم فرش فده 
البته در طی این روند خطای ناگریز (۵) به‌وجود 
خواهد آمد. در روش پیشنهادی این مقاله سعی 
می‌گردد با استفاده از کاهش تکنیکی ماتریس ضرایب 
تعداد معادلات (معادلات تکراری یا وابسته) کاهش 
یابد. این روش به‌لحاظ مبانی ریاضی دارای دو مزیت 
ها نیمه ۲۱ و٩‏ رف )یال یسوط یسدنه 


ریا زاس کسیر رسای نا اخبتال دود 
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ارتقای کارآیی شناسایی پارامترهای فیزیکی سازه با استفاده از... 


بردارهای موازی (و تقریباً موازی) از بدوضعی مسئله 
(تکینگی‌ها)» خطای حاصل از این چرخش (۵) کمتر 
شود. در نتیجه خطای کمتری در نتایج مجهولات )0( 
به‌وجود خواهد آمد. ۲ کاهش تعداد بردارهای 
موازی و شبه‌موازی موجب کاهش وابسته‌های عددی 
و خحطی ماتریس ضرایب و درنتبجه کاهش عدد 
وضعیت آن و پایداری بیشتر پاسخ در حضور نویز 


خواهد شد. 


حیل ۲ 


شکل ۱ نحوهٌ عملکرد پیش‌شرطساز در یک دستگاه دو معادله دو 


مجهول که با چرخش یکی از بردارهای شبه‌موازی همراه است 


لازم به توضیح است روند پیاده‌سازی الگوریتم 
کاهش تکینگی در ماتریس بلوکی 13 به‌صورتی بوده 
است که در ابتدا اين الگوریتم روی ماتریس شتاب ‏ 
لحاظ شده است و بعد از شناسایی ستون‌های وابسته از 
ماتریس شتاب حذف گردیده‌اند. در ادامه به‌صورت 
متناظر این ستون‌ها از ماتریس‌های سرعت و جابه‌جایی 
و هم‌چنین ماتریس نیروی ورودی(۳) نیز حذف 
گردیده‌اند. 
لم: اگر ۸ ماتریس مربعی (8«0) باشد و ماتریس ظ 
بالاشاف له اویش با اشفاده ازرروشن توش باشتاه 


آن‌گاه: 
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است اگر همه درایه‌های سطر ام ماتریس 13 صفر 


تال 


روش بالامثلئی‌سازی 

او ریس تقوح ح ان سح ظر وک معت ها نی 

(+(6 3۰ ,620 2۰ و6 :1) (در ایین‌جا ء یک 

عملگر مقدماتی (سطری مقدماتی) است. و 4 دو 
ادا رسب ره رهق اک ایکا 

می‌باشد) فقط عمل دوم و سوم استفاده گردد. 

(جابه‌جایی سطر و یا ستون در ماتریس ۸ نداشته 

باشیم). (۷ 6 ۸) 

0 زز۸ ( > 1 1) به‌عنوان عنصر محورگیری 
انتخاب شود و قبل از اجرای اعمال سطری مقدماتی 
در صورت امکان ناصفر گردد. 

۲ در روند اجرای الگوریتم هرکجا ز:۸۵ برابر صفر شد. 
محورگیری از بر وپز۸ ادامه یابد. (ه > 1+ ) 
اثبات: فرض کنید اشکالرهای غیرصفره (ع..» :2۸2 رین) 
وجود داشته باشد و سطرهای ۸ نیز به‌صورت :۸ ,۸ 
ان رصن اه اسف را عع ۱ 
آن‌ گاه اگر 0 < :۸ + 2,۵ +۰۰۰ + وهی( + ۵ب 
باشد. سطر [ ام ماتریس ۸ وابستهُ خطی (به سطرهای 

ماقبل خود) خواهد بود. 

به‌منظور استفاده از الگوریتم فوق در حل معادلة 

)٩(‏ دو نکن مهم بایستی مد نظر قرار گیرد. اول این‌که 

ماتریس ؟ ماتریسی ۳ < 10 (7 << ۳) است که برای 

استفاده از قضيه فوق باید تبدیل به زیرماتریس‌های 

۲ گر دد. دوم این که حالتی که در الگوریتم فوق 

ا تاه ری ار و رب 2 

برای وابستگی سطر ام ماتریس ۸ درنظر گرفته شد. 

در صورتی است که ماتریس ۸ دارای «قص رتبه» 

باشد. به‌عبارتی سطرها و يا ستون‌هایی در ماتریس ۸ 


هستند که به سایر سطر ها و ستون های ماتریس ۸ 


سال بیست و هشتم شمارة دو ۱۳۹۵ 


۳۹ 


واگ ظظرر قا را با هکت هافر شبات 
فیزیکی بدوضع که مبتنی بر جمع‌آوری اطلاعات 
میدانی هستند از نوع مسائل «بدحالت» می‌باشند. 
ار این ی تساه کی اریز شا هلا ار اریز 
ع < :۸ به‌جای رابطة قبل مقیاس بدوضعی خواهد 
بود. درواقع در این حالت به‌جای وابستگی خحطی 
بردارهای سطری و ستونی ۸۵ وابستگی عددی آنها مد 
نظر قرار هی کیرد؛ بنابراین به‌جای استفاده از عدد صفر 
در تشخیص وابستگی در الگوریتم ذکر شده بایستی 
یک عدد نزدیک به صفر اختیار گردد. 


پررسی نتایج شناسایی بر روی سازهٌ هشت طبقه 
نتایج شناسایی به‌کمک روش ارائه‌شده بر روی 
سازه‌های سه و هشت طبقه و تغییر بارگذاری در 
طبقات. با هم مقایسه شده‌اند. این نتایج نشان می‌دهد 
هر چه درجات آزادی سازه بیشتر باشد باتوجه به 
بو کر خواهد فد تاتیر استفاده از کتاهشن نکنیگ 

ماتریس در بهبود نتایج بیشتر خواهد بود. این موضوع 
با مالت‌های مار بتا از ابخاط محل امال :ناو بر روا 
دو سازهةٌ سه و هشت طبقه مورد بررسی و تجزیه و 
حالات فوق به‌منظور مقایسه با نتایج غفوری آشتیانی و 


قاسمی [31] آورده شده است. 


تشریح سازه, نحوة بارگذاری و مدل‌سازی سیستم. 
سازةٌ مورد نظر یکک قاب دوبعدی هشت طبقةٌ غیربرشی با 
میرایی نامتناسب می‌باشد. با فرض صلب بودن سقف‌هاه 
درجات آزادی سازه و جرم طبقات در تراز تیرها متم رکز 
شده اند و این درجات آزادی به‌صورت افقی مستقل 
تعریف گردیده‌اند. ماتریس سختی سازه با استفاده از 


مدل اجزای محدود و ماتریس میرایی سیستم. با 


نشريةُ مهندسی عمران فردوسی 


جایگزینی ماتریس سختی سازه با هندسة مشابه ولی 
اعضای متفاوت محاسبه شده است. خحصوصیات سازه در 
(شکل ۲) آورده شده است. 

با توجه به فرکانس‌های سازه به‌عنوان نیروی 
ورودی به سیستم. بار هارمونیک جارویی (با 
فرکانس‌های متغیر ۱ تا ۱۰ هرتز) در طبقه اول لحاظ 
تلم سر طیق ان پاسخ شا اوه مه دسنتت. | نله 
است. به‌جهت مدل‌سازی نویز نیز مطابق شکل (۲ ,2 
به‌شکل یک نویز سفید با توزیع نرمال و مقدار میانگین 
صفر که میزان آن براساس نسبت مقدار جذر مربعات 
سیگنال نویز به‌مقدار میانگین جذر مربعات پاسخ در 
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ارام کارا ای را فا فیریکین سار تم اوه از 


درجات آزادی پاسخ می‌باشد به پاسخ‌هاو نیروی 


ورودی صحیح افزوده شده است. 
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شکل ۲ سازه نمونه و خصوصیات آن 


نشرية مهندسی عمران فردوسی 


سال بیست و هشتم شمارهٌ دو ۱۳۹۵ 


مجید قاسمی- بابک تقوی-سعید عباس بندی 


نتایج شناسایی 

حساسیت شناسایی سیستم در دوحالت استفاده و عدم 
استفاده از الگوریتم کاهش تکینگی ماتریس در تعیین 
بان فیس با استفاده آن قبهبانت‌های بسرعت 
و جابه‌جایی به‌میزان نویز (۱ تا ۱۰ درصد) در پاسخ 
شتاب و نیروی ورودی در شکل (4) مقایسه شده‌اند. 
خطای نشان داده‌شده در این شکل. براساس میانگین 
خطای درایه‌های قطری ماتریس‌های تعیین‌شدهٌ سازه در 
۰ نمونه‌سازی به‌روش مونت کارلو محاسبه شده‌اند 
که نتایج نشان می‌دهد برتری نسبی از نظر مقدار خطا با 
وفاتش گنه اه کاهشی ککتکن ادوس کنو که 
این برتری با افزايش سطح نویز نمایان‌تر می‌گردد. 
بتابراین کاهش خطای "۵ در اثر پیش شرط‌ساز تابع 
جریمه که به‌عنوان یک ایده و هدف مطرح شده بود در 
عمل مشاهده می‌گردد. دلیل اين‌که چرا بهبود نتایج با 
افزایش سطح نویز بیشتر می‌شود نیز به خصوصیات 
شناسایی سیستم در حوزه فرکانس برمی‌گردد. در واقع 
شناسایی سیستم در حوزه فرکانس در یک بازة خاص 
(بازةُ مفید فرکانسی) صورت می‌گيرد. با افزایش سطح 
نویز یک سری داده نامفید به اين بازه وارد می‌شود که 
در ماتریس ضرایب به‌شکل بردار نمود پیدا می‌کند. این 
بردارها با بردارهای موجود در ماتریس ضرایب 
متشابهند و موجب ناپایداری پاسخ و افزايش خطای ۸۸ 
می‌گردند. درنتیجه با حذف آنها به‌کمک روش کاهش 
تکینگی ماتریس بهبود نتایج به‌دست خواهد آمد. لازم 
به ذکر است داده‌های تامفید ذکرشده در سطح نویز 
پایین قادر به ورود به یازهُ مفید فرکانسی نیستند. 
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اف ۱۱۳ 


شکل ۶ مقایسٌ نتایج شناسایی عناصر قطری ماتریس‌های 
خحصوصیات سیستم براساس تغییرات میزان نویز در پاسخ‌ها و 
نیروی در حالت بارگذاری هارمونیک جارویی ۱ تا ۱۰ هرتز بر 
روی طبقه اول از سازه هشت طبقه در دو حالت: (2) بدون استفاده 
از کاهش تکینگی ماتریس. (0) با استفاده از کاهش تکینگی 
ماتریس(18510)؛ (در سه شکل دیگر نیز مقايسهةٌ نظیربه‌نظیر 
ماتریس‌های خحصوصیات سیستم در دو روش نشان داده شده 


است) 


نشریهة مهندسی عمران فردوسی 


۳۲ 


هدف دیگری که در اين ارائه مورد توجه قرار 
گرفته است کاهش ناپایداری پاسخ می‌باشد. در یک 
مسئله پایدار افزايش نسبی سطح نویز و نتایج حاصل 
رابطهٌ منظمی با یکدیگر دارند. همان‌طور که در شکل 
(۶4) ملاحظه می‌گردد نحوه اف]زایش خطاهمای 
خحصوصیات تعیین شده در آثر افزایش میزان نویز در 
ال کته از کاشتن تکتکی تایرشن امتانه سا 
باشد رفتار منظم‌تری نسبت به حالت غیر آن دارد که 
ان ی نها هه کی کی ها مرن 
پایدارسازی پاسخ (کاهش اثر بدوضعی در پاسخ 
مسئله) رسید. 


بررسی آماری دقت و قابلیت درایه‌های 
ماتریس‌های مشخصة سیستم 
همان‌طور که در شکل (۶) مشاهده گردید نتایج 
شناسایی برای عناصر قطری ماتریس‌های مشخصه 
سازه در هر دو حالت استفاده و عدم استفاده از کاهش 
تکینگی ماتریس, تا سطح نویز 4۵ تقریبا برابر 
بوده‌است و از این میزان نویز به بعد به‌تدریج کاراین 
روش کاهش تکینگی ماتریس در کاهش خطای پاسخ 
نمایان می‌گردد. در این قسمت برای اطمینان از قابلیت 
روش پیشنهادی نتایج آماری ۱۰۰۰ بار شناسایی برای 
کليء درایه‌های ماتریس‌های مشخصهء سازه با 
نمونه‌سازی مونت کارلو مربوط به سطح نویز ۳ (که 
نتایج شناسایی نزدیک به هم دارند) بررسی گردیده و 
در جداول (۱) و (۲) آورده شده است. بدین منظور دو 
معیار آماری متوسط خطا و ضریب تغییرات برای 
بررسی کاهش خطا و پایدارسازی پاسخ برای کلية 
درایه‌های ماتریس مشخصات سیستم سازه‌ای در دو 
حالت استفاده و علدم استفاده از کاهش تکینگی 
ماتریس محاسبه و با هم مقایسه شده‌اند. نتایج نشان 


می‌دهد درصد خطای عناصر قطری ماتریس‌های جرم. 


نشریهُ مهندسی عمران فردوسی 


از کارا ای پراش فا فیریکین سار موه از 


میرایی و سختی در هر دو حالت تقریباً یکسان, و برای 

عناصر غیر قطری برتری نسبی با روشی است که از 

کاهش تکینگی ماتریس استفاده می‌کند. بنابراین تا 
سطح نویز ۸۳ استفاده از کاهش تکینگی ماتریس تأثیر 
چندانی در کاهش خطای شناسایی ندارد. هم‌چنین 

ضریب تغییرات برای کلیةٌ درایه‌های ماتریس‌های جرم. 

میرایی و سختی (عناصر قطری و غیرقطری) در حالت 

استفاده از کاهش تکینگی ماتریس نسبت به حالت عدم 
استفاده از آن نتایج مطلوب‌تری دارد. بنابراین از کاهش 

ضریب تغییرات به‌عنوان یک معیار پایداربودن پاسخ. 

می‌توان فهمید در سطح نویز کم (همانند ۳/) نیز 

که تلو رسفا فان تقاط مي باشل بان دیکی 

نتایج مشاهده‌شده می‌توان به موارد زیر اشاره نمود. 

۱) در حالت عدم استفاده از کاهش تکینگی ماتریس, 
میانگین خطای کليٌ درایه‌های ماتریس‌های جرم 
سختی و میرایی به‌ترتیب ۱/۷۸ ۰/۷۲ و 1۱/4۲ 
تک 

۲ در حالت استفاده از کاهش تکینگی میانگین حطای 
کلیة درایه‌های ماتریس‌های جرم. سختی و میرایی 
به‌ترتیب 0.۱/۸۰ ۰/۵۷ و ۸۱/۱۳ می‌باشد. 

۳ در حالت عدم استفاده از کاهش تکینگی» میانگین 
ضرایب تغییرات درایه‌های ماتریس‌های جرم 
سختی و میرایی به‌ترتیب ۱/۱۲ ۸۰/۳۳ و 1۱/۱۵ 

یرصان یوار ان مکی فانک 
ضرایب تغییرات درایه‌های ماتریس‌های جرم 
سختی و میرایی به‌ترتیب 0.۱/۹۶ 0/۲۸ و 1.۱/۰۲ 
بات 

۵ به‌طور کلی در هر دو حالت دقت و قابلیت اعتماد 
درایه‌های غیر قطری ماتریس‌های تعیین‌شده بیشتر 
از عناصر قطری هستند. 


سال بیست و هشتم شمارهٌ دو ۱۳۹۵ 


مجید قاسمی- بابک تفوی-سعید عیاس بندی ۳۳ 


7 به‌طور کلی در هر دو حالت قابلیت اعتماد ماتریس وی یت ین اک از اش 
سختی بیش از ماتریس‌های جرم و میرایی می‌باشد 
جدول ۱ خصوصیات تعبین‌شدة ماتریس‌های میرایی و سختی در حالت عدم استفاده از روش کاهش تکینگی ماتریس با ایجاد ۳/ نویز در 


پاسخ‌ها و نیروی ورودی برای حالت بارگذاری هارمونیک در یک گره از سازه (درصد خحطای مقیاس‌شده و ضریب تغییرات به‌ترتیب از چپ 


به راست در پرانتز نشان داده شده است) 
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سال بیست و هشتم شمارهةٌ دو. ۱۳۹۵ نشریهٌ مهندسی عمران فردوسی 


۳ ازا کارا سای پراش فا فریکین سار اوه از 


جدول ۲ خصوصیات تعیین شد؛ ماتریس‌های میرایی و سختی در حالت استفاده از کاهش تکینگی ماتریس با ایجاد ۸۳ نویز در پاسخ‌ها و 
نیروی ورودی برای حالت بارگذاری هارمونیک در یک گره از سازه (درصد خطای مقیاس‌شده و ضریب تغییرات به‌ترتیب از چپ به راست 


در پرانتز نشان داده شده است) 
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نتیحه گیری درجات آزادی سازه تحت یک تحریک اجباری و 


در این مقاله با استفاده از اندازه‌گیری پاسخ سازه در کاهش تکینگی ماتریس‌های حاصل از پاسخ‌هاء 


نشرية مهندسی عمران فردوسی سال بیست و هشتم شمارهٌ دو. ۱۳۹۵ 


مجید قاسمی- بابک تفوی-سعید عیاس بندی 


ماتریس‌های مشخصة سیستم خطی (شامل ماتریس‌های 
جرم. میرایبی و سختی) بااستفاده از حل معکوس 
معادلات حرکت در حوزه فرکانس در روند شناسایی 
سیستم به روش به‌روزرسانی مدل تعیین گردیده‌اند. اگر 
پاسخ‌های سازه بدون نویز باشند. ماتریس‌های مشخصه 
سیستم خطی به‌صورت دقیق تعیین می‌گردند. ولی در 
حالت پاسخ‌های نویزدار بهینه ترین مقادیر ماتریس‌های 
خحصوصیات سازه بااستفاده از حذف تکینگی 
ماتریس‌های پاسخ. روش حداقل مربعات و کمینه 
نمودن مجموع نیروهای ماندگار معادلات حرکت در 
همه درجات آزادی و همه نقاط فرکانسی منتخب 
تعیین می‌گردند. در روش مذکور از کاهش تکینگی 
قسمت حقیقی ماتریس تبدیل فوریه شتاب سازه» تحت 
یک تحریک اجباری در حل معکوس معادلات حرکت 
استفاده شده است. 

به‌منظور بهبود بخشیدن حساسیت روش به‌میزان 
نویز» فرم جبری مسئله مورد بازبینی و بررسی قرار 
گرفته است و الگوریتمی پيشنهاد گردیده است که قادر 


راجع 


۳۵ 


به کاهش عدد وضعیت ماتریس نگاشت به‌کمک حذف 
تکینگی‌های این ماتریس می‌باشد. به‌منظور بررسی 
کارایی روش پیشنهادی در تعیین ماتریس‌های 
خحصوصیات سیستم. قاب هشت‌طبفهُ غیربرشی انتخاب 
کده اس هن حالت کلین ان واده کط زوین کناهتشن 
تکینگی ماتریس. با کاهش عدد وضعیت ماتریس پاسخ 
دقت خصوصیات تعبین شده را بهبود می‌بخشد. این 
قابلیت با افزایش سطح نویز نمایان‌تر می‌شود و می‌توان 
به‌کمک آن شناسایی‌های با برداشت‌های پراختلال را 
به‌نحو مطلوب‌تر انجام داد. درپایان به‌منظور بررسی 
آماری نتایج. دو معیار پذیرفته‌شدهة میانگین خطا و 
ضریب تغییرات در یک سطح نویز نسبتاً کم» برای 
کلیةٌ درایه‌های ماتریس‌های مشخصه سیستم مورد 
تجزیه و تحلیل قرار گرفت و نشان داده شد استفاده از 
کاهش تکینگی ماتریس می‌تواند در پایدارسازی پاسخ 
و کاهش خطا در روند شناسایی سیستم عملکرد 
ناسین خاقته باشتا: 


۲ 501186۲۲۸۷۱6۲۲ روم‌تتاا تا صز من)عع]ع(۲ ععفصون ۵ فل‌مطاع۷ متصصهصرو نار رنه تاوع ۷ رشرل6ه1۷]01۲ 
2008(۰) ,0۲6 ۷ 

5 6 م6066 ماد رقصمتاهتام۴ آهناصمزع] 1 اممتاتوظ و۲ وصصعاطا0:ظ 1۳۳7۵۲96 ر.ت16۲0 ۷ ,15210۷ 
1998(۰) ,1۷۲601 

40 احتتاععصا م۶ کاعتنمو رحصعنگ ر"عصعاط۳۵ظ . 0عع0ظ۱۱۳-۳] عاعمعزنا هه اصعتزام نات صفط۳ ).۲۱21۲66۲ 
(1998) ,120610۳2ظظ ری ۷]۵)68۵ 0عتاج۸ 

,0۵0ص نایدا صتقتوی لاو فم‌تفصه تعلتمصنه اعمصاع ۵۶ مموع7 مامت0 ۲و۸ .5 رحصتن۴! 
۶ ۷66ص اطع روعع‌معنم5 صمتای‌تصصصصمت جع صمتاعصصمص تم نت۳2 90.1 ۲6۵0۲۲ [هع1عع 1[ 
1990(۰) ,60۳001087 ]1 

۷۰ ۳۱ .۲ 51۸۲ ,"ومع اف دزمح-عص0 102 معصمتنل‌جوم-۱از «انته‌متداه ۵1 عوععل ۸" ,.9.۲,مصن۴ 
(1991) ,645-653 .00 ,(4) 12 ,۷۵۱ .۸۵۵۵1 

,(26)6 ,۷۵۱ .0۳۱0 ,50 .له ر۵06۲ ام 0و یی عاقتتمع ر.فرتطعنن .روط ر.آرهانه 
2005(۰) ,1955-1988 2۳۰ 


0مصمتانل‌جم-۱11 «اتت‌متهاته تم معط صمتای‌گنمع۷ معتمصاط ر.ورتطفنن .روط را ردالعون .رهاط 
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۳3 از ام کارا سای پراش فا فیریکین سار اهاز 


2005(۰) ,(2081656 15) 269-287 0۵۰ و(15)3 ۷۵۱ ,دول 1۳۵۳۶۰ رکصعاوی تقعطن[ 

۲ ۱60اتصطاتای ,2006 ر028۵65 41 امه اصتمم-مصتاهم۱ متقسهعم رو رتطعنن رآ رهالهن ,نت۴ 
[عنصطعع 1 .مرصیا/کصمن)متاطانام/ع0.ع تتاحا تفن 3 مر :۳ وم ماطهانههه وا عصنتمعتم رطمتاقناطانام 
۶ ازوته۲نصنا عسطاصه! روع‌ممعنمه ممتای‌تصصصومت 0ص ومتاقصصمگص به والنمع ۳‏ ر05.12 ۲۱۵۵0۵۲۲ 
2005(۰) ,60۳001087 1[ 

حللز تاتتمطنصاه هم م0مصاعصه صمتای‌گنه تما رو رتطفنن) ری روط ر.آیهالعن ر.آیرهاط 
20080656(,)2005(۰ 3 ط) ,269-287 00۰ و(3) 15 ,۷۵۱ ,5۸ .1۲۵06 ر"قصعاورو تمعصنا 0عصمتات جوم 

,ما۷۵ او رتافد معصمنل‌جم- از جع ان فصمتاهناوم »زه۲طععاه تهعصنا نامز ی ریکارت1 1۷۲052۳06 
.(1995) ,499-313 .00 ,(4) 22 ,۲۷۵۱ 

ممناه»ازو0 دمص هه صمتعصعمه صمت‌صن-تهآنوصنه عط م) 0عاق1ع۲ و0مطاعوص 0صه معط .م0 رصه50 
1974(۰) ,798-824 00۰ ,(4) 11 ,۷۵۱ رفظ ۱۷۷۸۸۸۵۲۰ ۲۰ 5120/۰ ,0ص اف عطع ۵۶ کصمتاجنوء لوتععاصا 101 
۵ ۱۵ عممتای نموه طان؟ ممتاهمعصه تملامطم م۲۱ ۲۵۲ صطیتنمعام اماام2۲ظ. رب‌ارتاشنآر]۷رتانآ 
(1997) ,89-101 .00 ,64 ۷۵۱۰ ,.طاق۱۷ ۵۱۱۵۵ .ل ۱ رصع2۲00۱ 

514۸۷۲ اف ج ۵۶ مونع۲هت00نعو 0ج مناج تهلنام‌صنه معط عصتاماتملمی رصمطی #ههر.۳۲ ی روتامی) 
.(1965) ,205-224 00۰ ,2 ۷۵۱۰ .۸۵81 ۲ 

و0۵00 فص مامطععنام ها صمتاصامی. قمتعتاوو تققع. تممطنا ینمی هه ربق ظرتممحتعتاظ 
1965(۰) 

مص رمتاتامی ومتعیاوو تفع 4صه مومع ماه تقمانام‌طف ری رط‌عطتم ص2۳ ریت ,0ات6 
1970(۰) ,403-420 00۰ 14 ۷۵۱۰ ر.طاه]1۷ 

ویک112 ۷ "صمانوهممعع0 ماو تقانم‌صته عطه مممصی وم ماه 0۲۵0وصصد مه .1 رصقفطن 
.(1982) ,72-83 ۲۳0۰ و(8)1 ,۲۷۵1 

هاگ 220 .50 .[ 512 روهظ امجمههنلزط ۵۶ ماه تمانمدته عامتنهم م۷۷ رصفطی . هه .گر اعع[ 
(1990) ,873-912 0۳۰ ,(11)5 ,۷۵۱ ,باتا0) 

۵ اک زمعاه لو مهن مه همه متفه ملظ راعاتهط0ه ر.۷ ,۳۲۵۲۱8۵۲00 
94(۰) ,191-229 ۲0۰ و(67)2 ۷۵1۰ ۷۲06۱6۲۳۵16( 

هد م1 تفت عمط ۵۶ عصمناعنوء احهعاصا املع متام مه 0مصاعصر همه م۸ رحعوصه]۲ 
(1971) ,616-622 0۳۰ و(8)03 ,۷۵۱ .قح .مصتاا .7 ]کر ععناله۷ تقانام‌دزه 

1 ما کومتای‌ناموج ان عصعنووی تهعمنا 0عصمتانصمع-۱1 که ممتابلمی یمه مصا صرق رن لرطهته ۲۷ 
257-267,)1973(۰ 0۳۰ «(10)2 ۷۵1۰ .۸۵81 .۲عصت _ .ل 51۸۷۲ رکعاصامزم 2۵560 

۱۱۱-0 عاع6ون0 تععصن تم ولمطاعصه ته‌زمعصته مصومی ۵ ممتافصتصمیه آه‌تامهتص هر رب لرط278 ۲۷ 
1979(۰) ,100-111 ۲۳0۰ و(21)1 .۷۵۱ ,16۷716 ]91۸77 ,"عصع2۲۵0۱ 

۰ ,16۷۱6۷ ]51۸7 0۲۷۲۵۲۰ع-ب] عطا ۵۶ وصحعص رها فصعاطمتج 0عومه-[1۱ عاعونل 0۶ ولو افص رت ۰ظ رجه‌وضه۲۱ 
.(1992) ,561-80 .0۳0 ,(34)4 


۱۱۱-0 019066 12۵00801 ۲ه۱نعع عط رز ع۲رناعسیا عمط ۵۶ عونا مط1 .ظن1 تهعی۱) رین .ظ ۲127565 
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مجید قاسمی- بابک تفوی- سعید عیاس بندی 


1993(۰) ,1487-503 .00 ,14 ۷۵۱۰ عایام۵ .56 [ 51۸2 ."عصصع2۲۵۵۱ 

تهعصنا ۵۶ عصعنووه ع۲هنوو 0عمتانصم- ۱ ۵ متام ممتامتمم-صنظ ۱ رتم۷ ۲۷.2۵0 کارط(۵01۵0 ۲۷ 
2000(۰) ,119-124 0۳0۰ و(13)7 ۷۵۱۰ اعب ۷۲۵۰ ۸۵۵۱۰ ,"عصه‌تاهناوع 

۱۵00۰ صمتاه2تهناعع۲ مج همه فصماطمتص معاهلنمه ولمعتمعصز ۵ ممتانتمک ر.لافرتمصمصل 1 
1963(۰) ,1035-1038 .00 ,4 .1۳010 ,1۷۲205 ]50۷16 ,501-504 و(151)1963 991و رد2 .۸۵۱6۵0 
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2006(۰) ,5877-5890 

۲ طمتامناتافجمع۲ مصجاهصاصز موزع۲ظ۳] ری رتاحاحصقط رل رل‌وداعذاط ر.ظ ر۷]2265 ما ,1۷612 .۷ ,ظع1([)۷ 
۰ ,(13)5 ,۲۷۵ ,2005 م۲۵۵6 مه معممعنمک ص همه موع۲ه1 رد موی مه رتم2۵۳ 
2005(۰) ,507-543 

جچ هصنعدا ععتتمناتی ۳۳۵60 ۵۶ ومتای‌گتاجع0] مامه۳۲۵ ۷۷۰. رم۱۷۸۷۵۱۱۸ 2۳0 ر.ظ رتالع20170] و.ظ بلا۷ 
,619-639 00۰ ,36 1۳۰ )ناگ م۲8 علمنوط نهد ,۱۷۲6۵۸۵0 مصتاهل من-اع۱۷۵0 غصعصصعا۴ عاتصز۳ 11۳۳2۴۵۷۵۵ 
2007(۰) 

۳۵۳۵ ,۷/۵۵ ,۷۵۱۵۲ .ل ی فتاههونا [ع۱۷۲۵۵ طد صمتاهم2هانامم ,]۱۷۲۰ رالع۲ ۳۶ 20 ر.ظ رکتاتتااز1 
2008(۰) ,440-478 .00 ر(01:75)4 ۲۷ 

۶ صمنانام5 عوزع۲ظ1 هصنونا ود ۱۷۱6۵0 جمتایتاع10 وروی ,هقی هه بصعت طف ۳210۲1 
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۳۸ ارفا گرا این را فافر فزریکین شفره اوه زین 
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